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The conserving approximation scheme to many-body problems was developed by KadanofT and 
Baym using the functional-derivative approach. Another approach for the Hubbard model also 
satisfies conservation laws, but in addition it satisfies the Pauli principle and a number of sum rules. 
A concise formal derivation of that approach, using functional derivatives, is given in this conference 
paper to highlight formal analogies and differences with conserving approximations. 

I. INTRODUCTION 

Although quantum many-body theory is in principle exact, any practical application of the formahsm necessarily 
leads to the violation of some exact results. For example, it was noticed long ago that whenever infinite subsets of 
diagrams are resummed, conservation laws will be violated if certain rules are not followed. As another example of 
a possible inconsistency, one can obtain a relation between number of particles, chemical potential and temperature 
from either a Green function or from a derivative of the free energy. How can we make sure that the diagrammatic 
series for the Green function and that for the free energy give the same result? Both of these problems, namely 
obtaining response functions that satisfy conservation laws and making sure that the free energy is "thermodynamically 
consistent" with the Green function, can be solved using an approach that was developed by Kadanoff and Baym [1] 
[2] many years ago. This approach, described in the paper by Bickers (Chapter 6) in these Proceedings, allows one to 
obtain so-called "conserving approximations" [3] that have been extensively used in the last decade. However, as in 
any approximate calculation, some exact identities are not satisfied by conserving approximations. [4,6] For example, 
the Pauli principle is violated, integration over coupling constant of the potential energy obtained from susceptibilities 
does not lead back to the starting free energy, and the irreducible verticies generated by a functional derivative of the 
self-energy are not the same as the vertex corrections entering the self-energy itself. [7] In addition, since there is an 
infinite number of choice of diagrams for the Luttinger-Ward functional, there is also an infinite number of conserving 
approximations and it is difficult to decide on the validity of an approximation based on the criterion of conservation 
only. Further criticism of this approach has appeared in Refs. [4] and [5] 

In this paper, we wish to present a new derivation of a non- diagrammatic approach to the repulsive Hubbard model 
that was developed a few years ago. [4] An analogous derivation in the case of the attractive Hubbard model was 
presented at this workshop by Bumsoo Kyung and has already appeared. [8,9] Here, we go back to our earlier results on 
the repulsive Hubbard model, but we use a formal approach that highlights analogies and differences with conserving 
approximations. This paper will strive to be concise with the explicit purpose of being clear without getting lost into 
details. The published derivation for the attractive Hubbard model contains more of these details. [8] 

In section II we first present a derivation of conserving approximations that can be extended, in section III, to 
the derivation of our non-diagrammatic approach. This derivation focuses first on response functions and on the 
self-consistent determination of irreducible spin and charge vertices. That determination is done in such a way that 
the Pauli principle in its simplest form is satisfied. In other words, crossing symmetry (Chapter 6) is satisfied only 
for field operators on a single site, all at the same time. This is the so-called two-particle self-consistent (TPSC) 
part of the approach. We then improve our approximation for the self-energy by using the TPSC results in the exact 
non Hartree-Fock part of the self-energy. Vertices and Green function entering the final expression are at the same 
level of approximation. No Migdal theorem is assumed. The effects of transverse and longitudinal ffuctuations are 
considered in separate subsections before the final expression can be written down. The internal-accuracy check [4] is 
also presented and used in the derivation. We end up in Sec. IV with a discussion of extensions of the approach that 
are now being developed. The appendix addresses a technical detail. 
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II. FUNCTIONAL DERIVATIVE FORMALISM AND CONSERVING APPROXIMATIONS 



In this section, we first discuss single-particle properties, then response functions and, finally, we present the Hartree- 

Fock approximation as an example of a conserving approximation. The latter will bo helpful as a starting point for 
the derivation of the two-particle self-consistent (TPSC) approach in the next section. 



A. Single-particle properties 

Following functional methods of the Schwinger school [1,2,10], we begin with the generating function with source 
fields (pa and field destruction operators '0 in the grand canonical ensemble 

\nZ[4>] = In TV [e-^i"-^^)^, (e-^W)'^^?^)^'^))] (1) 

We adopt the convention that 1 stands for the position and imaginary time indices (ri , ri ) . The overbar means 

summation over every lattice site and integration over imaginary-time from to 0. T^ is the time-ordering operator. 
The propagator in the presence of the source field is obtained from functional differentiation 

Ga (1, 2; m = - {i,, (1) (2))^ = -^^y (2) 

Prom now on, the time-ordering operator in averages, (), is implicit. Physically relevant correlation functions are 
obtained for {(/)} = but it is extremely convenient to keep finite {(/)} in intermediate steps of the calculation. 

Using the equation of motion for the field ip and the definition of the self-energy, one obtains the Dyson equation 
in the presence of the source field [11] 

(Go 1 - G = 1 EG ; = 0(7^ - - S (3) 

where, from the commutator of the interacting part of the Hubbard Hamiltonian H, one obtains 

(1, T; {<!>}) (T, 2; {<j>}) = -U {i^l, (l+) (1) (1) i^l (2))^ • (4) 
The imaginary time in 1+ is infinitesimally larger than in 1. 



B. Response functions 



Response (four-point) functions for spin and charge excitations can be obtained from functional derivatives {dG/6(j)) 
of the source-dependent propagator. Following the standard approach and using matrix notation to abbreviate the 
summations and integrations we have, 

GG"^ = 1 (5) 



5G 



G-^ + G 



SG- 



= 0. 



(6) 



Using the Dyson equation (3) G ^ = Gq 



S this may be rewritten 



5G _ ^ SG- 
S(f> S(j) 



G 



G~G + G^G, 



(7) 



where the symbol - reminds us that the neighboring labels of the propagators have to be the same as those of the 
(j) in the functional derivative. If perturbation theory converges, we may write the self-energy as a functional of the 
propagator. From the chain rule, one then obtains an integral equation for the response function in the particle-hole 
channel that is the analog of the Bethe-Salpeter equation in the particle-particle channel 



dG 
1^ 



= G-G + G 



SESG 

sgJ^ 



G. 



(8) 
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The labels of the propagators in the last term are attached to the self energy, as in Eq.(7). [12] Vertices appropriate 
for spin and charge responses are given, respectively, by 

_ ^ _ ^ . jj _ ^ , ^ (n\ 



C. Hartree-Fock and RPA as an example 



As an example of calculation of response functions, consider the Hartree-Fock approximation which corresponds to 
factoring the four-point function in the definition of the self-energy Eq.(4) as if there were no correlations, 

(1,1; {4>]) Gf (T, 2; {</,}) = UG^l, (l, 1+; {cl>}) G^ (1, 2; {.^}) . (10) 
Multiplying the above equation by (G^) ^ , we are left with 

(1, 2; {</.}) = C/G?, (1, 1+; {</.}) (5(1-2), (11) 



<5Ef (1,2;{0}) 



^Gf (3, 4; {,/.}) 



f/5 (1 - 2) 5 (3 - 1) 5 (4 - 2) , 



(12) 



which, when substituted in the integral equation (8) for the response function, tells us that we have generated the 
random phase approximation (RPA) with, from Eq.(9), Usp = Uch = U. 



III. ANOTHER APPROACH 



The approach developed in Ref. [4], consists in two steps, corresponding to the following two subsections. First, one 
computes response functions from the TPSC approach, inspired by earlier work of Singwi. [13] Second, one generates 
an improved approximation for the self-energy starting from an exact expression for E that explicitly separates the 
infinite-frequency limit from the lower-frequency contribution. The method also contains an internal-accuracy check 
that is discussed in the last subsection. 



A. First step: two-particle self-consistency for G^^\'E^^\ rip = Usp and T^jj* = Uch 



In conserving approximations, the self-energy is obtained from a functional derivative S [G] =6^ [G] /(5G of $ the 
Luttinger-Ward functional, which is itself computed from a set of diagrams. To liberate ourselves from diagrams, we 
start instead from the exact expression for the self-energy, Eq.(4) and notice that when label 2 equals 1+, the right- 
hand side of this equation is equal to double-occupancy (n-fn^). Factoring as in Hartree-Fock amounts to assuming 
no correlations. Instead, we should insist that {n-ini) should be obtained self-consistently. After all, in the Hubbard 

model, there are only two local four point functions: (n-fn^) and (n^"^ = (j^f^ ■ The latter is given exactly, through 

the Pauli principle, by (ji'^'^ = (j^f^ — ('^t) — ('^i) — ^7^' when the filling n is known. In a way, (n-i-n^) in the 
self-energy equation (4), can be considered as an initial condition for the four point ftmction when one of the points, 
2, separates from all the others which are at 1. When that label 2 does not coincide with 1, it becomes more reasonable 
to factor d la Hartree-Fock. These physical ideas are implemented by postulating 

(1, T; {0}) G(i) (T, 2; {0}) = A{^|G« (l, 1+; {0}) G^^) (1, 2; {<P}) (13) 

where ^{<^} depends on external field and is chosen such that the exact result [14] 

S^. (1, T; {<^}) G, (T, 1+; {<!>}) = [/ (nt (1) (1))^ (14) 

is satisfied. It is easy to see that the solution is 
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- ''(nT(l)>,W(l)>/ 

Substituting A^^j back into our ansatz Eq.(13) we obtain the self-energy by right-multiplying by (gW) : 

S(i) (1, 2; {<^}) = ^{^}GW (1, 1+; {<^}) 5(1-2). (16) 

Wc arc now ready to obtain irreducible vertices using the prescription of the previous section, Eq.(9), namely through 
functional derivatives of S with respect to G. In the calculation of Usp, the functional derivative of (n-fn^) / ((n-f) (n^)) 
drops out, so we are left with, [15] 



{0}=o '^^T 



= ^,,j.o = C/7^. (17) 



The renormalization of this irreducible vertex may be physically understood as coming from Kanamori-Brueckner 
screening. [4] To close the system of equations, we need to know double-occupancy. It may be found self-consistently 
using the fluctuation-dissipation theorem and the Pauli principle. First notice that the Pauli principle, (n^) = (n^), 
implies that 

{n^-n{f) = {n^) + {n{)-2{n^n^) (18) 

while the fluctuation-dissipation theorem tells us that (^{n-\ — Jt-x)^^ is given by the equal-time equal-position 

imaginary-time susceptibility Xsp- Since x7p il) = Xo^il) ~ ^^sp with q = (q,27rnT), this is equivalent to the 
equation 

Xo{q) -n-2{n^n^) (19) 



-T — 



^UspXoiq) 



that, with Eq.(17) for Usp, gives us double-occupancy. 

The functional-derivative procedure generates an expression for the charge vertex Uch which involves the functional 
derivative of (ntn|) / ((^^t) i^l)) which contains six point functions that one does not really know how to evaluate. 
But, if we again assume that the vertex Uch is a constant, it is simply determined by the requirement that charge 
fluctuations also satisfy the fluctuation-dissipation theorem and the Pauli principle. 



N ^1 



-n + 2{npii) -n^. (20) 



Note that, in principle, S^^^ also depends on double-occupancy, but since S^^^ is a constant, it is absorbed in the 

definition of the chemical potential and we do not need to worry about it in this case. That is why the non-interacting 
irreducible susceptibility xo('7) appears in the expressions for the susceptibility, even though it should be evaluated 
with G(^) that contains S'^^. One can check that spin and charge conservation are satisfied by our susceptibilities. 

Detailed comparisons with quantum Monte Carlo simulations (QMC) [16,4,5,17] have shown that Eqs.(17)(19)(20) 
give predictions that are quantitative at the few percent level, in regions of parameter space where size effects are 
negligible. This remains true even at couplings of the order of the bandwidth and when second-neighbor hopping 
t' is present. [18] QMC size effects become important at half-filling below a crossover temperature Tx where the 
renormalized-classical regime appears. Even though equation (17) for Usp fails at n = 1, = in this regime, by 
assuming that (n|n|) is temperature independent below Tx one obtains a qualitatively correct description of the 
renormalized-classical regime. The universality class of our theory however is O (A^ = oo) instead of O (A^ = 3) . [19] 
A rough estimate of the renormalized chemical potential (or equivalently of S^^^), is given in appendix. 

B. Second step: an improved self-energy E^^^ 

Collective modes are less influenced by details of the single-particle properties than the other way around. We thus 
wish now to obtain an improved approximation for the self-energy that takes advantage of the fact that we have found 
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accurate approximations for the low-frequency spin and charge fluctuations. We begin from the general definition of 
the self-energy Eq.(4) obtained from Dyson's equation. The right-hand side of that equation can be obtained either 
from a functional derivative with respect to an external field that is diagonal in spin, as in our generating function 
Eq.(l), or by a functional derivative of (^ip-a (1) '/'J- (2))^ with respect to a transverse external field (pt- These two 
approaches will be considered in turn below. They give a self-energy formula that takes into account, respectively, 
longitudinal and transverse fluctuations. Crossing symmetry, rotational symmetry and sum rules will dictate the final 
formula for the improved self-energy E^^) that will be presented at the end of the subsection on the consistency check. 



1. Longitudinal channel 

The right-hand side of the general definition of the self-energy Eq.(4) may be written as 

SG. {1,2; {<!>}) 



S,(l,l) G„ (1,2) = -U 



{4>}=o 



(1,1+) G<.(1,2) 



(21) 



The last term is the Hartree-Fock contribution. It gives the exact result for the self-energy in the limit lo —>■ oo. [4] 
The 5Ga/54)-cr term is thus a contribution to lower frequencies and it comes from the spin and charge fluctuations. 
Right-multiplying the last equation by G~^ and replacing the lower energy part SGa/5<j)-a by its general expression 
in terms of irreducible vertices, Eq.(8) we find 



(1, 2) = UG^ll (1, 1+) S{l-2)- UG^^^ (1,3) 



5S« (3,2;{<ji}) 



(4, 5; {</.}) 



<5G« (4,5;{0}) 



{0}=O 



-.(1+,1) 



{<^}=o. 



(22) 



Every quantity appearing on the right-hand side of that equation has been taken from the TPSC results. This means 

in particular that the irreducible vertices ST,'^'^ /SG^^ are at the same level of approximation as the Green functions 

Gct ^ and self-energies In approaches that assume that Migdal's theorem applies to spin and charge fluctuations, 
one often sees renormalized Green functions G^^^ appearing on the right-hand side along with unrenormalized vertices, 

5E„/SG^, -> U. 

In terms of Usp and Uch in Fourier space, the above formula [5] reads. 
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(23) 



2. Transverse channel 



In the transverse channel, the calculation basically has to be redone from scratch. It is closely analogous to the 
attractive Hubbard model case, for which detailed calculations have been published. [8] We will thus be concise. The 
generating function in a transverse field is 

InZ = InTV [e-^i"-^^)^^ (e-^KT)^-(T.5)*i(2)-V.I(T)^+(T,2)v.T(5)^] (24) 
The corresponding spin-space matrix Green function 



G(l,2;{<^t}) = - ) / (25) 

' Vi(l)V' (2)y (^i(l)^ (2) ' 



obeys the matrix Dyson equation 

G-i (1, 2; {<^J) = Go-i (1 - 2) - S (1, 2; {4>t}) - (26) 

with 
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<j)- 



and 



S(1,1;{0J)G(1,2;{0J) = -C7 



JvI(l+)V'l(l)^T(l)^j(2) 
Jvj (l+)VT(l)V'i(l)V'| (2) 



V'|(l+)V'i(l)V'T(l)^I(2) 
V| (l+)VT(l)V'i(l)V'I (2) 



(27) 



(28) 



We need to obtain, in turn, the renormahzed vertex Usp Eq.(17), and finally find an improved formula for the 

sclf-cncrgy Starting from the matrix equation GG^^ = 1 and following the procedure below Eq.(5), one of the 

two response functions that do not vanish in zero external field obeys, 



5G2i(l,2;{0j) 



S<l>+ (3,4) 



= G..(l,3)Gn(4,2) + G..(l,2)^pHi(||^ 
{0j=o ^G2i (6, 7; {0t}) 



SG21 (ej-AM) 



t}=0 



S<l>+ (3,4) 



Gn (3,2) 



{4>t}=o 



(29) 



where the subscripts denote matrix elements in spin-space. There is an analogous equation for SG12 (1, 2) /S^- (3, 4) . 
Note that 

SG21 (1, 1) /(5</.+ (2, 2) = - {S+ (1) S. (2)) = - (V| (1) (1) i^l (2) Vt (2)) = -X+- (1, 2) • (30) 

The usual Hartree-Fock factorization for the self-energy would transform Eq.(29) into the RPA equation for transverse 
spin fluctuations {S+ (1) S- (2)). 

The value of S^^^ and of the corresponding vertices at the TPSC level are obtained using steps analogous to those 
in the longitudinal channel. We factor the self-energy in the Hartree-Fock manner. 



SW (l,T;{0t})G«(T,2;{</.J)=-A 



Gi'^{l,l+;{M) -G« (1, 1+; {</.*}) \ /g« (1,2; {.^t}) G« (1, 2; {</,,}) 
^^'^ ^ -G^V (1, 1+; {M) (1, 1+; {M) J \ G^V (1, 2; {4>t}) (1, 2; {M) 

(31) 

but we correct this factorization by the factor A^^^-y which is determined in such a way that when 2 — > 1+ the exact re- 
sult for the four point function on the right-hand side of Eq.(28) is recovered. Using (1^) ''Pi (1) V"? (1) V"} (1^)^ = 

(Vl (1+) V'l (1+) V'i (1) (1)) = and the analo gous result for the other off-diagonal element, the exact result for 
2^1+ is, 



E(^) (1,T; iM) G(^) (T, 1+; {M) = f/ ( (1) 



(32) 



Equating with the right-hand side of the approximate result Eq.(31) when 2 — > 1+ determines the value of A^^^y 
by a simple 2x2 matrix inversion. From this, one extracts the off-diagonal component of the self energy and the 
corresponding irreducible vertex that is needed for response functions in Eq.(29) 



5G^V(3,4;{0a) 



(33) 



Substituting in the equation for the response function Eq.(29) one precisely recovers for the transverse spin fluctuations 
the same TPSC result as for the longitudinal fluctuations: xjp (q) = Xo^{q) ~ ^Usp with Ugp given by Eq.(17). The 
determination of Usp through the fluctuation-dissipation theorem leads again to Eq.(19). We thus recover the results 
expected from rotational invariance. 

To move to the second level of approximation for the self-energy, we return to the exact definition for the self-energy 
Eq.(28) but this time in zero applied external field. The first diagonal component can be written as 



Ell (1,2) = [/ 



^G2i (l.2:{(/)a) 
6ct>+ (1+,1) 



GriM2,2). 



(34) 



t}=0 
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Using the exact result for the transverse response, Eq.(29) this takes the form 



En (1, 2) = UG22 (1, 1+) (5(1-2) + UG22 (1, 3) 



51221 (3,2;{.^,}) 



SG21 (6,7; {4.t}) 



SG21 (6,7;{<^J) 



{<^t}=0 



(35) 



i<t>t}=0 



where, as in the longitudinal case, the high-frequency Hartree-Fock result is now explicit. Substituting on the right- 
hand side the TPSC (i.e. level (1)) results, we obtain an improved approximation for the self-energy due to transverse 
spin fluctuations. In momentum space, it reads. 



(36) 



C. Internal accuracy check 



The equation (14) that relates Tr [EG] to double occupancy (potential energy) relates purely single-particle quanti- 
ties, E and G, to a quantity that may be computed from two-particle quantities, namely spin and charge correlation 
functions. In our approach, it can be checked [4] that ^Tr [E'^'G'^-'] is exactly equal to U {n^ni) , with {nini) 
computed at the TPSC level. More specifically, we find. 



iTr 
2 



E(2)g(i)1 = lim 



fe 



(37) 



One can use the difference between ^Tr [E^^^G*^^^] and ^Tr [E^^^G*^^-*] as an internal accuracy check of the theory. In 
the pscudogap regime at n = 1 for example, the breakdown of the theory is clearly indicated by the growing difference 
between the two quantities. Ref. [9] gives a detailed table that illustrates these facts in the case of the attractive 
Hubbard model. 

The Tr [Y.^^^G^^^] formula can also be used to help in the interpretation of the two different results obtained 
above for the self-energy, Eqs.(23) and (36). Without the approximation that the TPSC results should be used on 
the right-hand side, the results that follow from the corresponding exact expressions Eqs.(22) and (35) would be 
identical. To resolve this problem, we follow Ref. [22]. Fig. 1 of this paper shows the self-energy in terms of the 
fully reducible vertex T [q, k — k' , k + k' — q) . In both formulas for the self-energy, Eqs.(23) and (36), the dependence 
of r on the particle-particle channel center of mass momentum k + k' — q is neglected since this channel is not 
singular. The longitudinal version of the self-energy Eq.(23) takes good care of the singularity of F when its first 
argument q is near (tt, tt) . The transverse version does the same for the singular dependence near (tt, tt) of the second 
argument k — k' , which corresponds to the other particle-hole channel. One then expects that averaging the two 
possibilities gives a better approximation for F since it preserves crossing symmetry in the two particle-hole channels. 
Furthermore, one can verify that the longitudinal spin fluctuations in Eq.(23) contribute an amount U {n-^ni) /2 



to the consistency condition ^Tr ^5];^^^G(^^^ = U (n^nj^) and that each of the two transverse spin components also 

contribute U {n-[n\) /2 to |Tr ^E^^^^^G^^^^ = U {n-^n^) . Hence, averaging Eq.(23) and the expression in the transverse 

channel Eq.(36) also preserves rotational invariance. In addition, one verifies numerically that the exact sum rule [4] 
— / div' Im [E^ (k,w')] /tt = U'^n-„ (1 — n_o-) determining the high-frequency behavior is satisfied to a higher degree 
of accuracy. We thus obtain a self-energy formula that we called [22] "symmetric" 



= Un-a + 



U T 
'8N 



(38) 



l^cr {k)sym is different from so-called Berk-Schrieffer type expressions [20] that do not satisfy [4] the consistency 
condition between one- and two-particle properties, |Tr (EG) = U {n-[ni) . 



IV. DISCUSSION AND EXTENSIONS 



The approach described above is valid in the weak to intermediate coupling regime. It involves two steps. First, a 
self-consistent determination of double-occupancy and renormalized vertices that enter the dynamical susceptibilities 
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of the most important four-point correlation functions (density-density and spin-spin correlation functions). This is 
summarized by Eqs.(17)(19)(20). Conservation laws, such as charge and spin conservation, are satisfied. These results 
are then used to obtain an improved approximation for single-particle properties through the self-energy, Eq.(38). An 
internal accuracy check allows one to decide on the validity of the results in cases where QMC or other exact results 
are not available as references. The approach satisfies the Pauli principle, the Mermin- Wagner theorem, contains 
Kanamori-Brueckner screening and does not assume a Migdal theorem in the calculation of the self-energy. 

In addition to provide an accurate calculational tool for the Hubbard model, this methodology has allowed to 
develop insight into the Physics of this model. The physically most important result obtained to date is probably 
the detailed description of the Physics of the pseudogap that appears in t,wo dimensions in the renormalized-classical 
regime {Tiuj <C ksT) when the antiferromagnetic correlation length (superconducting correlation length in the case 
of the attractive model) becomes larger than the single particle thermal de Broglie wavelength. [5,4] This Physics 
explains a possible route to the destruction of the Fermi liquid in two dimensions. These results have been confirmed 
by extensive QMC calculations. [21,22,17] The pseudogap (depletion near w = 0) in Im (k^^-, w) appears along 
with precursors of the Bogoliubov quasiparticles (finite u peaks) of the ordered state. To extrapolate more deeply in 
the pseudogap regime, where strictly speaking the above method fails, one assumes that double-occupancy becomes 
temperature independent below the crossover temperature Tx where one enters the renormalized-classical regime. 

The above methodology has been applied successfully to the attractive Hubbard model. [8,9,23] It trivially applies 
to the Hubbard model with an arbitrary hopping matrix. S^^^^ can also be used to obtain consistent thermodynamic 
predictions. [24] Extensions to multiband problems are non-trivial but are being developed. [25] Recently, an extension 
of this approach has been used to obtain quantitative results for the spin and charge susceptibilities in the attractive 
Hubbard model. [26] Proceeding along the same lines for the repulsive Hubbard model, pairing correlations have 
been calculated. [27] They yield a dome shape dependence on doping of the d— wave superconducting transition 
temperatiirc Tc- The decrease of Tc near half-filling comes from the detrimental effect of opening a pseudogap. 
Phenomenological extensions to more complicated models with (i-wave superconductivity and antiferromagnetism 
have also been proposed. [28] Future directions include generalizations to the pseudogap regime, to states with broken 
symmetry, [29] to longer range interactions and to impurity models. It woiild also be extremely valuable to obtain 

the irreducible vertices 5T.)j' /dG]^/ that are consistent with the best estimate of the self-energy. 
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and Engineering Research Council (NSERC) of Canada and the Fonds pour la formation de Chercheurs et I'Aide a la 
Recherche (FCAR) of the Quebec government. We thank the Centre dc rccherches mathematiques for its hospitality. 
A.-M.S.T. holds a Tier I Canada Research Chair in Condensed Matter Physics. 



APPENDIX A: AN APPROXIMATE FORMULA FOR S^^' 



The equation Eq.(14) that was used to obtain finally Usp = U (n^ni) / ((n^) (n^)) would be different if we had taken 
the limit 2^1^ instead of 2 1+ in the general expression for SG, Eq.(4). More specifically, at {(f>} = 0, 



S,, (1, 1) G, (1, 1+) = U {ni (1) nt (1)) ; (l, l) (l, 1") = U (n^ (1) (n^ (1) - 1)) . 



(Al) 



Any approximation for the self-energy that has Hartree-Fock as its infinite-frequency limit will be such that the 
difference between the above two results, (l, l) G^- (l, 1"'") — S<7 (l, l) G^ (l, 1~) = U (n^ (1)) is satisfied. The 
proof, Eq.(44) of Ref. [8], is as follows: 



T 
N 



EE 7 



S (k, ikn) 



ikn - (£k - /u) - S (k, ik, 



U /g-ife„o- _ g-i/c„o+A 

n) ikn / ^ ' 



AT mm 



U{ni). 



(A2) 



In this expression, the first sum vanishes because we have added and subtracted a term that makes it convergent at 
infinity without the need for convergence factors e"*''"'^ . Hence, only the last sum survives. 
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Since (1,2) is a constant times (5(1 — 2), one can obtain two estimates of the constant depending on which of 
the two equations in (Al) one starts from. By analogy with Eq.(A2) one expects that the difference between these 
two estimates is related to the high frequency behavior of the true result while the average of the two estimates is 
related to the low-frequency behavior. This average is 

T) 1 

-{U + UsA^-n))j^. (A3) 

The exact result, U/2, is recovered at n = 1. For other fillings, the above formula gives a very rough estimate of the 
chemical potential shift induced by interactions. For example, for up to 4f on 8 x 8 lattices and temperatures of 

order t/4 in energy units, one finds results that deviate by up to 20% from the QMC results. The corresponding 
procedure in the attractive Hubbard model case [8] seems to work better. 
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